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Section  1 
INTRODUCTION 


The  Navy  relies  on  sonar-system  performance 
predictions  for  many  applications  —  from  concept  and  design 
through  direct  Fleet  tactical  aids.  Each  application 
involves  a  specific  amount  of  ensembling  over  conditions  and 
requires  a  certain  level  of  accuracy.  Performance  predic¬ 
tions  are  usually  generated  through  detection  or  tracking 
simulations,  employing  computer  models  for  engagements 
between  targets  and  sensors.  Whether  a  Monte  Carlo  or 
analytic  approach  is  used,  the  statistical  results  depend 
principally  on  the  scenario,  the  system  response  and  the 
acoustic  environment.  This  study  concentrates  on  the  last 
item,  the  environmental  acoustic  (EVA)  parameters  and  their 
fluctuations,  and  how  it  impacts  performance-prediction 
capability. 
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A  substantial  amount  of  Navy  R&D  effort  has 
been  concerned  with  understanding  and  modeling  the  averaged 
and  detailed  properties  of  EVA  parameters  (specifically, 
sound  transmission,  ambient  noise,  and  reverberation).  The 
critical  questions  related  to  this  work  are: 


(a)  Under  given  ensembling  conditions  (applica¬ 
tion),  how  sensitive  are  performance  model 
predictions  to  the  EVA  inputs? 

(b)  What  is  the  level  of  accuracy  (detail)  needed 
in  the  estimation  of  EVA  data  as  a  function 
of  the  type  of  ensembling? 
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(c)  Given  answers  to  (a)  and  (b),  how  should 
the  EVA  properties  and  system  performance  be 
estimated? 

“(d)  In  particular,  are  there  efficient  means  for 
calculating  measures  of  effectiveness, 
detection  statistics,  and  detailed  properties 
of  the  EVA  parameters? 

The  four  questions  are  non-trivial  and  cannot  be  resolved 
completely  in  this  investigation  (see,  e.g.,  Reference  1-1 
for  a  sample  of  recent  research  and  evidence  that  the 
acoustic  community  is  not  ready  to  address  the  problem) . 
The  present  effort  does,  however,  consider  important  parts 
of  each  question.  The  work  has  proceeded  in  two  phases,  as 
described  next. 

The  first  phase  (Ref.  1-2)  concentrated  on  low- 
frequency,  passive  sonar  and  questions  (c)  and  (d);  that  the 
Navy  had  developed  and  was  using  performance  prediction 
models  based  on  a  specific  approach:  random-process  simu¬ 
lations  for  acoustic  or  system  fluctuations  and  Monte-Carlo 
calculations  of  MOE ' s • (detection  probabilities,  holding 
times) .  It  was  determined  that  the  accuracy  of  the  random- 
process  approach  was  strongly  dependent  on  the  use  of 
accurate  input  parameters: 

•  Mean  Transmission  Loss  at  Array  Output  (TL). 

•  Mean  Noise  at  Array  Output  (Beam  Noise). 

•  Distributions  and  Temporal-Correlation  Proper¬ 
ties  of  TL  and  Beam-Noise  Fluctuations. 


It  was  further  concluded  that  the  random-process  approach 
could  be  improved  by  (i)  the  selection  of  random-process 
models  with  properties  more  like  those  of  actual  acoustic 
fluctuations,  and  by  (ii)  more  efficient  means  for  deter¬ 
mining  the  EVA  statistical  inputs,  such  as  those  listed 
above. 

The  second  and  current  phase  of  investigation 
has  two  concerns: 

I.  Initial  consideration  of  questions  (a) 
and  (b),  especially  as  related  to  calculating 
a  specific  MOE:  cumulative  detection  proba¬ 
bility  <CP/D>. 

II.  Improvement  of  the  random-process  approach 
in  response  to  the  recommendations  just  men¬ 
tioned  (i  &  ii). 

For  item  I,  the  first  priority  was  to  gain  an 
understanding  of  the  sensitivity  of  performance  predictions 
to  the  accuracy  of  EVA  inputs  for  various  levels  of  ensembl- 
ing  (over  target  tracks,  diurnal  variations,  etc.).  A 
special  case  in  which  a  signal-to-noise-ratio  (SNR)  detector 
operates  against  deterministic  SNR  was  used  to  study  the 
ensembling  question  in  great  detail.  Deterministic  and 
various  types  of  random  tracks  were  considered.  The  results 
showed  sensitivity  of  discrete  and  cumulative  detection 
probabilities  (Pq  and  CPd)  to  ensembling,  and  were  then 
extrapolated  to  more  general  cases.  Although  the  solution 
is  not  complete  (this  was  an  initial  and  limited  effort), 


the  illustrations  of  Section  2  should  be  of  considerable 
value  to  performance  modelers  for  determining  input-accuracy 
requirements  and  for  interpreting  their  predictions. 

The  second  problem  pursued  under  item  I  is  that  of 
practically  estimating  CPq  for  "realistic*'  acoustic 
environments.  Most  detection  simulations  (see  again  Refer¬ 
ence  1-20  employ  random-process  models  for  signal  and  noise 
fluctuation  with  special  forms  (e.g.,  exponential  autocor¬ 
relations)  which  generally  do  not  yield  faithful  facsimiles 
of  the  usual  quasi-periodic  variations.  Moreover,  there 
seems  to  be  some  confusion  about  the  underlying  properties 
of  signals  and  noise,  and  how  they  impact  instantaneous 
detections  (e.g.  via  ROC  curves),  and  how  they  affect 
overall  system  performance  (e.g.  CPj)).  Section  3  docu¬ 
ments  an  initial  effort  to  address  these  problems,  first  by 
carefully  describing  the  detection  problem  in  terms  of  basic 
signal  and  noise  properties,  and  then  by  deriving  an 
approach  for  calculating  CPp  for  a  "realistic"  fluctuating 
signal-to-noise  ratio  (SNR).  Again,  the  results  should  be 
valuable  to  performance  models,  both  as  aids  to  interpreting 
their  results  and  as  improvements  to  modeling  methods. 

Turn  now  to  the  second  item  (II)  of  interest  here: 
improvement  of  the  currently  popular  random-process  approach 
for  fluctuation  modeling  in  performance  analyses.  The 
effort  (just  described)  to  use  realistic  SNR  for  CPp 
calculations  represents  a  contribution  in  Problem  Area  II. 
However,  additional  work  concentrated  on  approaches  for 
obtaining  the  key  features  of  TL  and  noise  fluctuations 
without  the  extensive  computer  computations  usually  required. 
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Section  4  and  a  comparison  report  attack  the  TL 
fluctuation  problem.  The  premise  is  that  fluctuations  in 
transmission  for  most  passive  sonar  applications  are  caused 
by  relative  source/receiver  motion. 

A  separate  document  (Ref  1-3)  provides  a  solution 
to  a  problem  of  great  interest  in  the  EVA  community — how  to 
efficiently  model  TL  for  a  moving  source.  The  conclusion  is 
that,  under  reasonable  conditions,  variations  in  multipath 
interference  induced  by  source/receiver  motion  can  be 
accurately  modeled  with  ray  or  wave  theory  in  a  "frozen" 
ocean.  It  is.  not  necessary  to  integrate  the  time-dependent 
wave  equation.  The  resulting  shortcut  is  of  considerable 
importance,  and  the  analysis  is  summarized  in  this  document. 

With  the  cause  of  the  TL  fluctuations  and  an  ap¬ 
proach  to  model  them  identified,  an  efficient  ray-theoretic 
method  has  been  developed  to  provide  the  important  spatial 
properties  of  the  TL  fluctuations,  which  in  turn  are  con¬ 
verted  to  temporal  properties  for  dynamic  scenarios.  The 
technique  required  a  very  modest  amount  of  ray-tracing  and 
has  been  applied  to  both  bottom  and  waterborne  paths.  The 
computer  program  and  associated  graphics  package  have  been 
used  in  Section  4  to  test  the  results  against  detailed  TL 
model  predictions  and  to  investigate  the  sensitivity  of  the 
fluctuation  statistics  to  the  environmental  inputs.  It  is 
expected  that  the  algorithm  will  find  extensive  applications 
in  providing  inputs  to  performance  models. 

Just  as  for  TL,  a  short-cut  for  predicting  key 
features  of  beam-noise  fluctuations  at  low  frequency  (ship¬ 
generated  noise)  would  be  a  valuable  asset  for  performance- 


prediction  applications.  The  various  noise  fluctuation 
schemes  in  current  use  were  reviewed  in  detail  (from  Refer¬ 
ence  1-4).  Limitations  in  the  most  promising  candidates 
were  identified,  and  a  variation  of  an  existing  technique 
selected. *•  Results  were  tested  against  "Brute  Force"  simula¬ 
tions.  The  investigation  is  documented  in  Section  5. 

Finally,  given  TL  and  noise  fluctuation  proper¬ 
ties,  the  remaining  problem  under  item  II  is  to  identify 
random-process  models  which  more  closely  reflect  those 
properties  than  currently  popular  processes  (e.g. ,  Jump, 
Gauss-Markov,  Ehrenfest).  Section  6  recommends  certain 
auto-regressive  schemes  as  candidate  fluctuation  processes. 

In  summary,  this  study  has  considered  a  variety 
of  questions  related  to  sonar  performance  modeling.  Pro¬ 
gress  has  been  made  in  describing  the  general  sensitivity 
and  approach  for  calculating  CPj)  as  well  as  in  developing 
methods  for  generating  EVA  inputs  to  detection  and  tracking 
simulations. 
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Section  2 


DETECTION  PROBABILITIES  FOR  A 
REALISTIC  ACOUSTIC  ENVIRONMENT 


2.1  Introduction 

Sonar  system  performance  is  traditionally  gauged 
in  terms  of  instantaneous  detection  and  false-alarm  proba¬ 
bilities  (PD  and  Pfa)*  Among  the  most  important  mea¬ 
sures  of  effectiveness  for  both  surveillance  and  tactical 
ASW  systems  is  the  long-term,  accumulated  effects  of  numer¬ 
ous  detection  opportunities:  the  cumulative  detection 
probability  (CPp). 

Efforts  to  understand  and  model  CPj)  in  the  past 
have  been  principally  focused  on  two  approaches:  the 
independent  glimpse  model  and  random-process  simulations. 

In  the  first,  the  detection  opportunities  are  assumed  to 
become  statistically  independent  after  some  waiting  time,  so 
that  CPd  is  a  simple  function  of  individual  Pd's.  The 
second  is  concerned  with  the  acoustic  signal  and  noise, 
treating  each  (or  their  ratio)  as  structured  random  pro¬ 
cesses  (e.g.,  Markov,  with  Gaussian  marginal  densities). 
Detections  then  occur  whenever  thresholds  are  exceeded.  In 
both  cases  it  has  been  standard  procedure  to  model  the 
acoustic  signal  and  noise  as  very  smooth  functions  (of  time 
or  range)  to  which  are  added  random  fluctuations. 

Seldom  are  the  underlying  structures  of  the  signal 
and  noise  and  their  relationship  with  the  instantaneous  and 
cumulative  detection  process  properly  defined.  Concepts  of 


independence,  decorrelation,  and  Pd  itself  are  often 
confused,  especially  when  there  is  ensembling  over  target  or 
other  conditions.  Subsection  2.2  proposes  a  unified  view  of 
the  acoustic  detection  problem,  including  the  detector  and 
fluctuations. 

Note  next  that  research  in  environmental  acous¬ 
tics  over  the  past  several  years  has  led  to  a  new  under¬ 
standing  of  the  temporal  and  spatial  variations  of  both 
signals  and  ambient  noise.  In  particular,  it  is  usually  the 
case  that  the  fluctuations  are  quasi-periodic ,  with  predict¬ 
able  periods  and  amplitude  statistics,  and  that  there  are 
statistically  well-defined  states  of  slowly  varying,  smooth, 
mean  values.  Part  of  the  present  investigation  endeavors  to 
develop  algorithms  for  efficiently  predicting  these  important 
features  of  the  acoustic  environment  (see  Sections  4  and 
5).  This  section  deals  with  the  problem  of  calculating 
CPd  from  such  a  description  -  extending  the  "glimpse"  and 
"random-process"  approaches  to  the  more  complex  and  struc¬ 
tured  signal  and  noise  fields. 

Subsection  2.3  defines  the  general  CPq  problem, 
given  the  underlying  structure.  Traditional  characteriza¬ 
tions  of  the  signal  and  noise,  and  approaches  to  CPd 
calculations  are  summarized  in  Subsection  2.4.  Finally,  a 
"realistic"  model  of  the  acoustic  variables  and  a  means  to 
obtain  the  attendant  CPd  are  presented  in  the  final 
section . 
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2.2  Underlying  Structure  for  Detection  Process 

A  review  of  the  standard  references  for  acoustic 
detection  suggests  two  levels  of  detail.  From  a  signal- 
processor's  point  of  view  (e.g.  Reference  2-1  or  2-2), 
estimates  of  signal-plus-noise  and  noise  are  somehow  mani¬ 
pulated  (filtered,  integrated)  and  compared,  with  detection 
called  according  to  a  preselected  rule  (threshold).  The 
input  variables  are  modeled  as  random  processes  and  both 
correct  detections  and  false  alarms  are  possible  (with 
probabilities  Pj)  and  Pfa)*  It  is  standard  practice  to 
treat  one  detection  opportunity  (integration  period)  at  a 
time,  and  to  ignore  the  concept  of  CPp.  On  the  other 
hand,  sonar  performance  and  engagement  modelers  (see,  e.g., 
References  2-3  and  2-4)  tend  to'  assume  a  binary  detector 
operating  on  the  random  process  SE(t),  usually  the  inte¬ 
grated  signal  excess  estimated  from  the  signal-processor’s 
output.  A  correct  detection  occurs  when  SE(t)>0  and  does 
not  otherwise.  Pp  is  then  the  probability  that  SE>0,  and 
the  false  alarm  rate  is  not  explicitly  given.  Cumulative 
probability  (CPq)  is  then  the  chance  that  SE(t)>0  for  some 
t  in  an  extended  time  interval. 

The  two  views  are  not  very  consistent,  and, 
although  cases  can  be  fabricated  to  bring  them  into  agree¬ 
ment,  serious  complications  arise  when  the  signal-processor's 
detector  and  the  engagement-modeler's  random  process  must  be 
incorporated  in  a  careful  calculation  of  CPd.  How  does 
decorrelation  time  fit  in?  What  are  independent,  the 
detection  events  or  the  random-process  samples? 
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To  properly  treat  the  CPd  problem  a  scheme  for 
viewing  the  detection  process  is  proposed  next,  with  the 
intent  of  including  the  two  viewpoints  discussed  above. 

First,  let  X  and  Y  be  vector-valued  variables 
which  represent,  respectively,  the  received  signal  and  noise 
fields,  and  the  detection  rule  for  a  predetermined  detection- 
opportunity  interval.  Then,  in  general,  detection  proba¬ 
bility  (Pd)  is  the  probability  that  X  and  Y  take  on 
certain  values;  Pp^  has  a  similar  definition.  This 
almost  trivial  structure  is  intended  to  include  nearly  any 
signal  processor  scheme. 

The  complicating  and  crucial  second  step  is  to 
force  the  dependence  of  Pd  and  Pp^  on  X  and  Y  in  a  special 
form; 

Assume  the  existence  of  a  summary  description 
(statistic),  SN,  for  X  and  Y  at  each  time  interval 
such  that  Pd  and  PpA  depend  only  on  that 
description.  In  other  words  detections  and  false 
alarms  occur  according  to  values  assumed  by  SN  at 
each  opportunity. 

This  is  the  usual  approach  for  construction  of  an  ROC  curve, 
with  SN  some  averaged  signal-to-noise  ratio. 

Next,  if  SN  is  deterministic  (as  is  usually  the  case  for 
signal  processor's  ROC  curves),  then  Pp  and  Pfa  are  real 
numbers.  If,  however,  SN  is  a  random  variable,  then  so  are 
Pd  and  PpAJ  this  is  the  point  at  which  confusion  can 
begin.  Return  to  the  original  definition  of  Pd,  where  the 
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implication  is  that  the  probability  is  based  on  the  popula¬ 
tion  X  and  Y.  It  follows  that  if  SN  must  be  modeled  as  a 
random  variable  to  include  the  effect  of  those  populations, 
then  the  correct  value  of  Pq  will  be  obtained  only  as  an 
ensemble  over  the  population  of  SN.  Such  ensembling  takes 
the  form  of  an  expected  value  over  the  possible  values  (Z) 
of  SN: 

Pd({SN})  =  <Pd(SN|SN=z)>2  (2-1) 

This  type  of  expression  is  generally  true  for  any  random 
variable,  e.g.,  in  terms  of  the  mean  value,  W,  of  variable 

W 

P(W>C) 

=  P(W>C|W=x)  P(W=x)  dx 

=  E^(P( W>C | W) )  (2~2 

Continue  now  with  the  concept  of  cumulative  detection 
probability  (CPd(T) :  the  probability  that  over  a  fixed 
time  interval  (say  IO,T],  including  a  number  of  detection- 
opportunity  periods)  there  is  at  least  one  detection.  For 
the  first  definition  of  Pq,  it  is  the  probability  of  X  and 
Y  assuming  certain  values  over  at  least  one  detection 
period.  Likewise  for  SN^  or  SN(  *),  the  values  of  SN  at 
the  i-th  opportunity  interval 


o 


If  SN^  are  deterministic,  then  so  are  the 
associated  detection  probabilities,  Ppi. 
Moreover,  it  is  usual  to  presume  that 
depends  only  on  SN^,  and  not  on  any  previous 
or  future  values  of  SN.  In  that  case, 


CP  (T)  =  1  -  II  d-PDi)  (2-3) 

i 

the  usual  form  of  CPq  for  "independent" 
detection  opportunities. 

o  If  SN^  are  random  variables,  then  so  are  the 
Poi.  If  the  variables  are  independent, 

then  I 

CPD(T)  1  -  n  (1-E(PD.)) 

(i.e.,  we  have  ensembled  over  all  values  of  (SNj, 

SN2 

o  If  the  SNi  are  random,  but  not  necessarily 

Independent,  then  CPp(T)  is  non-trivial,  and  no 
general  formula  can  be  given  without  additional 
assumptions  about  SN^  . 

Subsequent  subsections  are  concerned  with  the  latter  problem. 

2.3  The  CPd  Problem 


The  aim  here  is  to  investigate  cumulative  detection 
probability  (CPd)  defined  as 
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m 


CP[>(T)  =  Probability  {Detection  Occurs  at  Least 

Once  in  the  Time  Travel 
[0 ,  T] } 


It  will  be  convenient  to  use  the  complement: 


Q(T)  =  1  -  CPq(T)  =  Probability  |No  Detection  in  [0,T]}. 


To  make  the  notation  somewhat  more  precise,  define 
the  following  events. 


A(t)  =  Event  that  detection  occurs  at  time  t, 

A(t,t)  *  Event  that  detection  occurs  at  least  once 
in  time  interval  (t,  t  +  t  ), 

A[t ,  t  ]  - - [t,  t  +  t  ], 

A(t)  and  A( t , t ) :  events  complementary  to  A(t)  and 

A(t ,  t  )  . 


Then 

PD(t)  =  P(A( t ) ) , 

CPD(t)  =  P( A(0, t ] ) , 

Q(t)  =  P(A(0,tJ). 
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It  is  common  practice  to  calculate  CPq  stepwise 
from  the  general  recursive  relationship: 


Q(T  +  At)  =  Q(T)  •  (1  -  P(T ; At ) ) , 


(.2-4) 


CPd(T  +  At)  •  1  -  (1  -  CPd(T))  (1  -  P(T ; At  )  ) , 


(2-4') 


=  CPD(T)  +  P(T ; At )  (1  -  CPD(T)) 


where 


P(T;At)  =  Probability  that  detection  occurs  in  the  interval 
'  (T,  T  +  At],  given  that  no  detection  has  occurred 
in  [0,T] 


=  P( A(T, AT ]  I  A[0,T] ) . 


(2-5) 


Equation  (2-4)  is  then  no  more  than  the  definition  of  the  condi¬ 
tional  probability.  When  Q(t)  =  0,  P(t;At)  is  undefined,  but 
may  be  arbitrarily  set  to  0  or  1  or  whatever.  If  the  underlying 
process  is  Markov,  then  A[0,TJ  can  be  replaced  by  A(_t)  for  any  t 
in  [0 ,T] . 


If  we  consider  discrete  time  steps  0,  At,  2At,  ..., 
NAt  =  T,  then 

Qn+1  =  Q( (N+l ) At )  =  Q(NAt)  •  (1  -  P(NAt ; At ) ) 


=  QN  *  (1  -  P(NAt ; At ) ) 


N  v 

...  =  /  77(1  -  P( JAt ; At )))  •  Qi  (2-( 


where  1-Q^is  the  "initial"  probability  of  detection  over 
(0,  At]. 


If  the  process  is  treated  as  continuous  in  t,  then 
the  existence  of  the  limit 


lim 

At+0 


^P(t;At)^ 


P(t) 


(2-7) 


implies  that 


Q'(t) 


lim  (Q(t+At)  -  Q(t ) 
\  At 


=  -Q(t)  •  p(t). 


(2-8) 


Then 


T 

Q(T)  =  Q(0)  •  exp  (-  f  p(t)dt), 

0 


(2-9) 


or 


T  „ 

CPD(T)  =  1  -  (1  -  CPD(0))  exp  (-^  p(t)dt).  (2-9*) 


In  some  cases  (2-9)  is  preferred  to  (2-6);  however  (2-7)  may 
not  exist  and  the  integration  time,  At,  is  usually  finite. 
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Whichever  applies,  discrete  or  continuous  time 
steps,  the  problem  is  to  calculate  P(t;At)  of  (2-5).  The 
remainder  of  this  section  is  concerned  with  models  of  SN(t), 
the  detection  process,  and  the  calculation  of  CP[)  using 
(2-5). 


2.4  Examples  and  Traditional  Approaches 


This  section  considers  progressively  more  diffi¬ 
cult  problems,  with  popular  models  of  SN(t)  and  the  detec¬ 
tion  process.  It  leads  toward  a  realistic  model  for  the 
case  of  interest,  to  be  attacked  in  the  next  section. 


2.4.1  Problem  1 


SN(t)  is  a  deterministic  function  of  t. 


Pj)( t )  =  P(SN(t)^L)  for  some  fixed  L. 


pFa ( t )  =  0 


Then 


PD(t)  = 


CPD(T)  = 


0  if  SN( t )  <  L 
1  if  SN(  t )  >,  L 

0  if  SN(t)  never  exceeds  L  for  0  <  t  <  T 
1  otherwise 


max  { Pp ( t ) }  . 
0  1  t  <  T 
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Note  that  Pq  and  CPd  are  completely  deterministic.  There 
is  no  "accumulation  of  independent  opportunities,"  and  CPd  is 
just  the  largest  Pd. 

Example  Problem  1  is  so  simple  that  the  use  of  (2-9)  or 
(2-6)  is  not  required. 

2.4.2  Problem  2 


SN(  t )  is  a  deterministic  function  of  t. 

Po(t)  =  G(SN(t)  i'*  v  deterministic  function,  but 
not  limited  to  values  {0,1}. 

PFA(t)  is  similarly  defined  [say,  H(SN(t))]. 


Although  it  appears  only  superficially  more  general 
than  Problem  1,  this  case  encompasses  many  of  the  difficulties 
associated  with  defining  and  calculating  CPd»  It  is  commonly 
encountered,  with  the  most  typical  example  being  that  in  which 
the  detector  performance  is  defined  by  an  ROC  function,  so  that 
Pq  depends  on  short-term,  averaged,  output  statistics  of 
signal-plus-noise  and  noise  (called  "SN(t)M  here). 

It  must  be  emphasized  at  the  outset  that  SN(t)  is  a 
deterministic  function.  Hence,  all  uncertainty  in  the  detec¬ 
tion  process  is  reflected  in  the  function  G  which  maps  the  real 
number  SN  into  a  detection  probability.  To  this  point,  there 
is  no  concern  about  any  underlying  signal  or  noise  or  process¬ 
ing  or  threshold;  it  is  all  contained  in  G  (and  the  Pp^ 
function).  However,  as  we  proceed  to  the  calculation  of  CPq» 
notice  that  the  lack  of  a  better  description  of  SN  and  Pd 
causes  difficulties. 
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At  each  point  or  interval  in  time  (discrete  or 
continuous)  there  is  an  opportunity  to  detect  quantified  by  the 
instantaneous  and  deterministic  P[>(t).  This  is  not  enough 
information,  however,  to  calculate  CPq.  In  particular,  the 
dependence  of  the  detection  events  from  one  time  step  to  the 
next  is  required. 


In  terms  of  detection  events,  we  know 
PpCt)  .  P( A( t ) )  or  PD(t)  =  P( A( t , At } ) , 

CPD(T)  =  P(A[o,T ] ) . 

If  the  events  A[o,At],  A(At,  2At ],..., A(NAt , (N+l ) At ] 
were  mutually  independent,  then 


P(kAt ; At )  =  P(A(kAt,At) |A[o,kAt]) 
=  P(  A(kAt ,  At ) ) 


and  from  (2-6), 


CPD  (N+l) At) 


1-P(  jAt  ;At) 


)) 


( 1-P(A ( o, At J 


=  1 


!N 

n 

j= i 


(1-P(A( jAt, Atl) 


( 1-P(A( o , At ] ) ) 


(2-10) 
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This  is  the  "independent  glimpse"  result,  so  often  applied  in 
surveillance.  If  on  the  other  hand,  all  of  the  cited  events 
were  totally  dependent,  then 

P(k At ; At )  .  P(A(kAt, At ]| A[o,kAt ]) 

=  0,  for  all  k, 

and  again  from  (2-6), 


CPD(  (N+l )  At  )  =  p(a[  o,  At  ])  =  p(A(  jAt ,  At  ]  ) ,  for  all  j. 

(2-11) 


For  a  continuum  of  time  steps,  the  events  A(t)  are  of 
interest.  For  the  case  of  independence  among  all  events,  the 
existence  of  p  from  (2-4), 

/  P(A(t ; At ) )  \ 

p(t)  *  -lim  1  .  /  , 

At-*-o  \  x  t 

gives  the  result  analogous  to  (2-10),  from  (2-9'): 

T 

CPd(T)  =  l-(l-CPD(0))exp(  -  (t)dt).  (2-12) 

A 

For  totally  dependent  events,  p  =  0,  and 


CPd(T)  =  CPD(0). 


(2-13) 
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These  special  cases  should  illustrate  the  importance  of 
defining  the  degree  of  dependence  among  the  events  A.  This 
goes  far  beyond  the  specification  of  instantaneous  Pp’s. 
Finally,  note  that  Pp  and  CPp  are  real  numbers,  not  random 
variables.  We  hear  too  often  in  the  context  of  this  problem 
the  meaningless  statement  that  the  Pq's  are  independent  or 
dependent . 


2.4.3  Probl em  3 


!SN(t)  is  a  random  process 

Pp(t)  =  P(SN(t)^L),  for  fixed  threshold  L 
PpA  is  fixed,  independent  of  SN 


This  is  perhaps  the  model  most  often  used  in  ASW 
detection  problems.  It  is  like  Problem  1  in  that  the  detec¬ 
tor  is  easily  defined  as  a  thresholding  process;  however,  the 
quantity  tested  is  no  longer  deterministic,  but  rather  a 
random  variable  at  each  time  step.  In  Problem  2,  the 
uncertainty  was  concentrated  in  the  detection  process;  here 
it  is  contained  in  the  underlying  input  to  the  detector  (the 
detector  behaves  deterministically).  The  level-crossing 
properties  of  the  random  process  SN(t)  are  what  determine 
Pq  and  CPq. 

As  an  aside,  we  note  that  confusion  can  arise 
when  the  deterministic  SN(t)  of  Problem  1  is  complicated,  and 
is  thus  replaced  with  a  simpler  random  process  which  is  sup¬ 
posed  to  incorporate  the  important  statistical  features  of  SN . 
The  trouble  is  that  such  descriptors  as  decorrelation  time  or 
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one-dimensional  moments  which  fully  determine  the  random 
process  may  not  be  appropriate  for  the  deterministic  time 
series  (e.g.,  a  nearly  sinusoidal  signal). 

Substantial  efforts  have  been  devoted  over  the  past 
ten  years  to  the  calculation  of  CPp  (and  other  MOE's)  for  a 
class  of  random  processes,  especially  Gauss-Markov,  Jump,  and 
related  variations  (Ref.  2-5  to  2-7).  When  the  random  process 
represents  the  fluctuation  term,  it  is  added  to  a  deterministic 
term  (e.g.,  the  smooth  signal),  and  the  level-crossing  CPp 
problem  is  non-trivial.  A  recursive  solution  for  the  Jump 
process  is  available  (Ref.  2-7),  but  no  such  general  result 
is  known  for  the  Gauss-Markov  or  others.  Since  analytic 
solutions  are  rare,  it  is  common  practice  to  employ  Monte 
Carlo  simulations  in  the  search  for  CPq. 

Example  To  illustrate  the  complexity  of  this  problem, 
consider  first  the  simplest  "constant-threshold”  case 
for  the  Gauss-Jump  process.  Let  SN(t)=X(t)  be  a  zero- 
mean  Jump  (A, a)  process  (Ref.  2-7),  and  let  detection 
occur  when  SN(t)  ^  o.  Then 

CPD(T)  =  P(X(t )>o  for  some  te  [o,T]) 

=  1  -  P(X(t )<o  for  all  t) 

=  1-/2  P(X(t)<o|k  jumps)»P(k  jumps) 

Ik 

1  -  (l-Po)e-ATPo,  (2-14) 

where  Po  =  P(X(t)>o)  for  fixed  t.  In  terms  of  (2-5) 

P  (t;  At)  =  i_e-AAtPo,  (2-15) 
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and 


lim  P(t;At)  =  p(t)  =  APo . 


(2-16) 


At-*-0  At 

Hence  (2-9^)  yields  (2-14). 

However,  if  SN(t)  is  modeled  as 
SN(t)  =  L(t)  +X(t) 

where  X(t)  is  as  above  and  L(t)  is  deterministic,  the  problem 
becomes  substantially  more  difficult.  For  the  simplest  case 
of  L(t)  monotone  decreasing,  and  with 

P(t)  =  PD(t), 

p(t)  =  XP(t)  -  (P’(t)/(1-P(t). 

Hence,  CP[)(t)  =  1  -  (1-P(t))exp  f  P(x)dT^. 

When  L(t)  is  oscillatory,  endpoint  conditions  must 
with  a  recursive  scheme. 


(2-17) 

(2-18) 

be  matched 
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Problem  4 


SN(t)  is  a  random  process 


PFA(t) 


as  in  Problem  2, 
and  deterministic 


This  problem  includes  all  of  the  complications  of 
Problems  2  and  3:  uncertainties  in  detectability  are  the  result 
of  random  acoustic  parameters  (embodied  in  SN(t))  and  a 
random  dectection  process  (represented  by  G).  It  is  the  most 
realistic  model  of  those  presented  so  far  since  the  detector  is 
in  practice  (and  in  most  performance  estimates)  not  a  determi¬ 
nistic  thresholder  and  since  the  acoustic  signal  and  noise  are 
best  seen  as  random  fluctuations  added  to  deterministic  func¬ 
tions. 


Careful  interpretation  of  Pd  and  CPd  is  required. 

G  is  meant  to  map  the  random  variable  SN(t)  into  a  single  Pd; 
it  is  not  intended  that  Pp  be  random.  Hence  G(SN(t))  could  be 
P(SN(t)^L)  as  in  Problem  3,  or  some  other  ensemble  probability. 
Likewise,  CPp  should  be  an  ensemble  function  of  the  random 
process  SN(t)  and  not  a  random  process  itself. 

Mathematical  approaches  to  calculating  CPp  must 
employ  the  underlying  detection-event  structure  introduced  in 
Problem  2.  Thus,  as  in  that  case,  more  information  about 
dependence  than  is  stated  above  is  required.  This  might 
made  clearer  if  Pd  is  viewed  as  follows  (see  Subsection  2.2): 


/' 


PD^ t )  =  J  P (Detect ion  | SN( t )  =  x) • p(SN( t )=x)  dx 
E jp( Detect ion | SN(t) ) j • 


(2-19) 


Similar  expressions  for  P(t;At)  and  CPp(t)  suggest  that 
these  probabilities  can  be  treated  as  expectations,  over  the 
random  process  (SN(t)),  of  conditional  detection  probabilities. 
The  conditionals  in  turn  can  be  known  only  when  the  underlying 
dependence  among  detection  events  is  specified. 
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It  should  suffice  to  say,  without  examples,  that 
Problem  4  is  difficult  to  solve  analytically.  It  does  fore¬ 
shadow  the  realistic  model  presented  in  the  next  section. 

Table  2-1  displays  the  salient  features  of  each  of 
the  four  problems  treated  here. 


2.5  CPp  for  a  Realistic  Acoustic/Detector  Model 


2.5.1  Acoustic  Variables 


We  concentrate  on  passive  sonar  applications.  A 
rather  general  formulation  for  most  cases  of  interest 
is  as  follows: 

o  Detector  performance  at  discrete  time  step  t^ 
depends  on  an  integrated  variable,  SN(ti), 
at  the  output  of  the  be°mformer/processor  system. 

o  There  are  a  number  of  States  or  conditions  which 
strongly  influence  the  behavior  of  SN.  The 
chances  of  a  state  obtaining  and  the  time  in 
state,  etc.  are  described  by  a  stochastic 
process  or  by  a  deterministic  specification. 

o  Whenever  a  single  state  applies  for  a  number  of 
time  steps,  SN(t)  has  significant  structure. 

This  structure  is  well-described  by  periodicity 
(or  quasi-periodicity)  of  sinusoid-like  time 
series  with  a  distribution  of  peak  amplitudes. 


I 


Problem  SN(t)  Pd(0  Pfa(*)  CPo(t) 


1  Deterministic  P(SN(t)>L)  0  max  Pp(t) 

t 


2 


Deterministic 


G(SN(t)), 

H(SN(t ) ) 

Depends  on 

Deterministic 

dependence  of 

"detection  events" 

3  Random  Process  P(SN(t)>L)  Fixed 


Found  from  level¬ 
crossing  proper¬ 
ties  of  SN(t) 


4 


Random  Process 

G( SN( t  )  ) > 

H(SN( t ) ) 

Depends  on  both 

Deterministic 

dependence  of 

"detection  events" 

and  properties  of 

random  process 
SN(  t ) 


TABLE  2-1 


l 


I 
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All  of  this  is  suggested  by  the  structure  of  acoustic  signal-to- 
noise  ratio  in  which  the  principal  fluctuations  are  exhibits  of 
regular  multipath  interference  patterns  and  the  states  are  such 
conditions  as:  CZ  propagation,  BB  propagation,  beam-free  noise, 
beam-occupied  noise,  and  combinations  thereof.  Sections  4  and 
5  of  this  report  discuss  the  signal  and  noise  properties  in  some 
detail . 

2.5.2  Special  Problems 

As  in  Subsection  2.4,  we  will  attack  the  CPD 
problem  in  steps  of  increasing  difficulty,  beginning  with  a 
deterministic  SN(t)  and  threshold  detector  and  finishing 
with  the  most  general  case.  As  the  results  are  obtained, 
the  effect  of  added  complexities  and  randomness  will  be 
studied  -  so  that  efficient  algorithms  can  be  developed  at 
appropriate  levels  of  detail. 

Problem  A  /SN(t)  deterministic^ 

jpD  =  P(SN(t)  £  L). 

lpFA  =  0 


Just  as  for  Problem  1,  CPpCT)  =  max  {PpCt)}, 

0^tv<T 

and  hence  either  0  or  1  depending  on  whether  SN(t)  ever 
exceeds  L.  When  SN(t)  is  a  quasi-sinusoid  in  each  state, 
we  need  only  find  the  highest  "peak"  in  each  state,  and 
then  the  greatest  of  these.  The  problem  is  trivial,  and 
in  fact  solutions  are  given  in  the  next  section  in  the  study 
of  sensitivity. 
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Problem  B 


SN( t )  deterministic,  quasi-sinusoid 
|  for  each  of  several  states 

jpD  =  G(SN),  deterministic 

PFA  =  H(SN) 


Just  as  for  Problem  2,  some  statement  must  be  made 
about  the  dependence  of  detection  events.  Suppose  first  that 
they  are  all  independent,  then  correlated.  This  problem  is 
nearly  equivalent  then  to  A,  and  the  approach  the  same  as 
that  described  in  Section  3. 


Problem  C 


SN(t)  a  random  process  within  each] 
\  state;  states  deterministic  J 


jPD  =  P(SN  >,  L)  ( 

|pFA  fixed  I 

This  is  a  difficult  problem,  but  a  solution  has 
been  found  for  a  realistic  special  case.  SN(t)  is  presumed 
to  consist  of  a  deterministic  sinusoid  modulated  by  a  jump 
process.  The  recursive  scheme  of  Ref.  2-7  has  been  applied, 
and  the  resulting  algorithm  tested. 
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Section  3 


SENSITIVITY  OF  PERFORMANCE  PREDICTIONS 
TO  ENSEMBLIN6 

In  the  previous  section,  great  care  was  taken  to  define  the  underlying 
signal  and  noise  processes  and  their  implications  for  Pr,  and  CPQ  calculations. 
Even  for  the  case  of  a  stable  source,  single  source/receiver  geometry, 
and  fixed  system  response  the  analysis  can  be  quite  complex,  requiring 
ensemble  averaging  over  fluctuating  parameters  and  non-trivial  accumulation 
of  instantaneous  detection  probabilities.  Nonetheless,  this  "deterministic" 
scenario  can  be  important  for  some  applications  (e.g.,  reconstruction  of 
a  controlled  experiment  or  explanation  of  a  detection  event  with  known 
target),  and  the  sensitivity  of  the  performance  model  outputs  to  the  details 
of  the  EVA-parameter  inputs  is  relatively  well  understood.  Turn  then  to 
other  common  apoli cations  in  which  the  single,  "deterministic"  event  is 
only  one  replication  of  many. 

It  is  the  rule  in  engagement  modeling,  and  often  the  case  in 
performance  analyses,  that  predictions  or  simulations  are  generated  for 
a  number  of  conditions  and  displayed  as  a  distribution  function  (or  its 
moments).  For  example,  a  simulation  model  might  yield  a  histogram  of 
detection  probabilities  determined  from  a  large  number  of  target  tracks, 
source  levels,  and  environmental  conditions.  The  obvious  questions  for 
this  study  are: 

•  Under  given  ensembling  conditions,  how  sensitive  are 

the  performance  predictions  to  the  environmental  acoustic 
(EVA)  parameters? 

Or,  in  other  words, 

•  What  is  the  level  of  accuracy  (detail)  needed  in  the 
estimate  of  EVA  data  as  a  function  of  the  type  (extent) 
of  ensembling? 


lb. 


This  question  is  addressed  below.  The  results  are  not  surprising  to  any¬ 
one  who  has  thought  about  the  problem,  but  have  implications  which  should 
be  (and  are  often  not)  taken  into  account  in  every  performance/engagement 
analysis. 

The  simplest,  most  elementary,  problem  is  treated  first  in  great 
detail.  The  conclusions  are  then  extended  to  more  general  cases. 

3. 1  Special  Case 

Assume  the  simplest  type  of  detector,  operating  on  a  quantity  to 
be  called  SN,  which  is  some  processed  function  of  the  signal  and  noise. 

This  quantity  SN  is  described  deterministically  as  a  function  of  range  (R) 
from  source  to  receiver,  and  detection  occurs  (for  certain)  when  SN(R)>0. 
The  source/receiver  geometry  as  a  function  of  time  is  the  ensembling 
characteristic,  and  it  is  embodied  in  the  range  variable  R(t).  Because 
of  the  simple  detector,  the  only  property  of  SN(R)  which  matters  is  whether 
or  not  it  exceeds  zero.  It.  can  thus  be  reduced  to  a  function  which  has 
values  either  0  or  1.  Suppose  that  SN(R)  has  only  finitely  many  changes 

/N 

in  value  over  the  ranges  0  to  R,  and  that  these  changes  occur  at  ranges: 

Rq  =  °,  Rlt  R2,...,  Rn_-,,  Rn  »  R. 

The  function  SN  can  then  be  defined  as: 


SN(R)  = 


0  for  Re  f  R^.  _  ^ ) ,  i  odd 
1  for  Re  IR^  ,  R. ),  i  even 


Finally,  let  L.  denote  the  length  of  the  intervals  { ,R^ ) 
SN(R) 


♦  LN 


rn-i  rn 
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Note  then  that  instantaneous  and  cumulative  detection  probabilities  are 
calculated  according  to: 

PD(t)  =  1  if  R(t)«L^  for  i  even 

CPD( t )  =  1  if  R(s)  «L.j  for  i  even,  for  some  s«  |0,t)  . 

If  R(t)  is  a  random  process,  then  so  are  PQ(t)  and  CPD(t),  and  we  con¬ 
sider  their  distributions  for  each  t  as  functions  of  the  distribution 
of  R(t) . 

Finally,  a  subtle  point  must  be  made.  Pg(t)  for  a  deterministic 
R(t)  is  simply 

PD(t)  =  P(SN(R(t) )  >  0) 

|  0  if  R(t)eL.. ,  i  odd, 

\  1  otherwise. 

If  R ( t )  is  a  random  variable,  then  we  define  Pp(t)  as  the  random  variable 

!0  when  R(t)£L. ,  i  odd 
1 

1  otherwise. 

which  takes  on  values  of  0  and  1  only,  with  a  distribution  induced  by  that 
of  R(t)  for  each  t.  What  some  might  call  Pp(t)  (viz.,  P(R(t)«L^,  i  even), 
is  in  our  case  the  expected  value  of  PD(t)  over  R(t), 

E(PQ(t))  =  l-P(R(t)*L.,  i  even) 

+  0*P(R(t)«L.,  i  odd). 

The  distribution  of  PD(t),  over  R(t),  is  given  by 

(  P(PD(t)  =  1 )  =  P(R(t)«Lr  i  even), 

(  P(P„(t)  «  0)  =  l-P(PD(t)  «  1). 

Similar  considerations  were  made  in  Section  2. 
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3.1.1  Random  Initial  Range,  Deterministic  and  Constant  Soeed 


The  relative  source/receiver  range  is  of  form 


R(t)  =  rQ  +  vt, 

where  v  is  constant  and  rQ  is  uniform  on  [0,  R^l.  To  eliminate  boundary 
effects,  we  suppose  that  V>0  and  SN(R)  is  continued  periodically  to 
R  = 

Calculation  of  the  basic  MOE's  is  then  straightforward. 


!1  if (r  +  vt)£L . ,  i  even 

o 

0  otherwise 


P  jPD(t)  =  1}  =  P|rQ  +  vteL. ,  i  even}  = 

/  N 

■  (&,•  /  £l,)  s  le/l 

i  even'  i  =  1 

E  {pdU)}=  lk*(LE/L)  +  0k(L0/L) 

=  le/l 

Var  (PD(t))  =  Le/L-(Le/L)2  =  (L^/L) * (LQ/L) 

(  1  if  r0eL^ ,  i  even 

CPD(t)  =  <  1  if  rQEL^ ,  i  odd,  and  (rQ  +  tv)>R^ 
y  0  if  r0el..,  i  odd,  and  (rQ  +  tv)  < R^ 
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P  { CPg(t)  =  O}  =  P(rQeL^ ,  i  odd;  and  (rQ  +  tv)*^-) 

;  5/  <Ri-i-<ro<Rf  - tv) 

:?.££)  ■ 


i  odd 

=  53  max 
i  odd 


(3-1) 


i  odd 
L.,-  >  tv 


A 

Let  {  L.j  |  be  a  reordering  of  the  odd-indexed  L.'s  such  that  for  odd 


Lj1Lj+2* 

Then  P|CPg(t)  =  0}  is  a  piecewise  linear  function  of  t 
with  nodes 


*  i./v,  j  odd;  P(CPD(t)  =  0) 


=  53 


L.-tv 


i  odd 
i>j 


and  slopes  of  -kv/L,  where  k  is  the  number  of  indices  with  L^>tv,  as 
shown  in  Figure  3-1. 

Of  course,  P|CPg{t)  =  l|  =  1-P  }CPg(t)  =  0},  which  is  again 

piecewise  linear  with  the  obvious  nodes.  Notice  that 

P  {CPD(0)  -  1}  -  le/l  »  P  (PD(t)  =11, 


P  {CPD( t)  *  1}  *  1  if  t>(L./v)  for  all  odd  i. 


If  Lz  =  L3  =  . 


for  all  j,  then 


L2j+1 


m(^) 


when  tv<Lj 


P  (CPD(t)  =  0}  * 


otherwise. 


Also,  in  any  case. 


mJo  [rlML 

I  L  ) 


<P  (CPD(t) 


=  0} < max 


0, 


m(wt) 


as  shown  in  Figure  3-1  as  dashed  lines. 


The  moments  of  CPQ(t)  follow  easily: 


E{(CPD(t) )k>  =  lk-  P(CPD(t)  =  1)  +  0k-P(CPD(t)  =  0) 


-  l-P(CPD(t)  *  0). 


3.1.2  Random  Speed,  Deterministic  Initial  Range 
As  above,  let 


R(t)  =  rQ  +  vt.  . 

However,  now  fix  rQ  and  let  v  be  uniformly  distributed  on  the  interval 
v-j.Vgj,  where  0<v-|<v2*  It  is  clear  that  pD(t)»  P{Pg(t)=l),  and  hence 
1/ 

E{Pd  (t)}  are  the  same  as  before. 
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Likewise, 


CPD(t) 


1 ,  if  rQ  (  Li ,  i  even 

1,  if  r0€l-j*  i  odd>  and  r0+tv  >  Ri 

0,  if  rQ  eL.. ,  i  odd,  and  rQ+tv  <  R_.. 


However, 


P{CP0(t)=0} 


0,  if  r^L^  i  even 


),  if  rQfL^,  i  odd.  and  ro+tv.j  >  R^ 

1),  if  r0tLi »  1  °dd. 

ro  +  tv2 


1,  if  rQ€L.. ,  i  odd,  and  rQ+tv2  <  R^ 


(3-2) 


Thus  P{CPD(t)=0}  as  a  function  of  t  is  either  the  constant  function 

zero  (if  r  <L,  i  even),  or  it  consists  of  three  seqments  -  the  constant 

01R-r  Rr  Rr 

function  1  for  t  <  J _ o  ,  hyperbolic  for  i  o  <  t  <  i  o  ,  and  the  constant 


nM 

function  zero  for  t>Jj _ o  .  This  is  shown  in  Figure  3-2. 

vi 
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3.1.3 


Suppose  now  that  rQ  and  v  are  each  uniformly  distributed. 
Once  again  PQ( t) ,  P{PD(t)=1},  E{P0k(t)},  and  CPQ(t)  are  the  same  as 
before.  Recalling  the  result  for  a  fixed  speed,  (3-1),  we  have 


P{CPD(t)  =  0} 


j  P{CPD(t)=0|  v} 


P  (v)  dv 


1 

v2-vl 


A  ^ 

Li”tv 


dv 


1 


'2  V1 


L^<tV 


Li-1st(v1+v2) 

L 


i 


E 

Li<tv 


As 


L 


dv 


(3-3) 


To  carry  this  further,  consider  three  cases: 

A 

Case  I:  0$t£  ^i.  ( i  • e  •  >  t.>tv,  for  all  i  and  for  all  v  € fvi  *v2  ) 

v2 

Since  there  are  no  t^tv,  (3-3)  becomes 
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P{CPD(t)=0}  = 


yi  Li-Ht(v1+v2) 

L^Htlv^Vg)  [ 


=  P{CPD(t)=0  !  v=is(v1+v2)} 


(3-4a) 


_  |/w  / 

This,  as  seen  in  section  3.1.1,  is  linear  with  slope  '|_. 


.  /v 

Case  II:  t  >  Hn  (i.e.,  L.ctv,  for  all  i  and  for  all  v*  v.  ,v, 

m  *  '  ^ 


Since  all  L^<tv,  (3-3)  becomes 


P{CPD(t)=0}  =  0. 


(3-4b) 


Case  III :  L1  <  t<  Lm 


For  a  given  t,  choose  j  such  that  L-  s„max  {L . } . 


Then  from  (3-3) 


J  V„ 


P{CPD(t) 


-01  -  22  V*<VV  -  TT^  •  £  f  ~ 

1  I  '2  1  1=1  7  (Li/t,v. 

^  max  '  1 


L.-tv 

i 


dv 


=  22  i-i~,'at(V1tV2l  .  22  - 


i=1  V2  - 


*-i  L- 

-max{r>vl}  .  Li-»at(v2+nBx{r,v1}) 

L 


Next,  choose  k  such  that 


Lk  *  min  {L1 } 

h>tvi 


Then 


P{CPD(t)=0}  = 


Li-Jst(v1+v2) 

L 


-52 


.  Lj  -  »»t(v2*L1/t) 
L 


( 3-4c) 


Finally,  it  can  be  shown  that  P{CP0(t)=0}  and  its  first 
derivative  are  continuous  for  t>0.  Higher  order  derivatives  in  general, 
however,  are  not  continuous.  Figure  3-3  shows  a  particular  P{CPg(t)=0}. 
The  vertical  dashed  lines  delineate  the  intervals  over  which  all  deriva¬ 
tives  are  continuous. 
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3.2 


Extrapolations 


The  special  case  given  above  was  examined  in  detail  for 
several  reasons: 

o  It  is  as  simple  a  case  as  can  be  contrived,  and  yields 
closed-form  solutions  (equations  (3-1),  (3-2)  and  (3-4)). 

o  The  form  of  the  solutions  allows  for  direct  examination 
of  the  sensitivity  of  the  detection  probabilities  to  the 
input  scenario  for  important  types  of  ensembling. 

o  With  some  modifications,  the  results  can  be  easily  extended 
to  more  general,  realistic  problems. 

Consider  next  some  observations  based  on  these  results. 

Note  first  that  SN(t)  is  determined  completely  by  the  interval 
lengths  |  }  and  the  range  rate  R(t).  When  R(t)  is  a  random  process, 

so  is  SN(t).  If  R(t)  is  Markov,  then  SN(t)  is  also  Markov.  Similarly, 
randomness  in  |  |  induces  randomness  in  SN,  for  deterministic  R(t). 

Next  consider  the  association  of  the  forms  of  R(t)  with  properties 
of  SN,  usually  identified  as  fluctuations: 


Special  Case 
Random  initial  range 

Random,  but  constant  speed 

Random  range  rate  ( R( t ) ) 


Extrapolation 

Imprecise  location  of 
peaks  and  nulls  in  SN(R) 

Imprecise  scaling  of  periods 
of  SN 

Random  periods  and  amplitudes 
for  SN 


Periodic  speed  Periodic  variations  in  SN 

Formulas  such  as  (3-1)  or  (3-4)  can  then  be  used  to  determine  precision 
needed  in  the  definition  of  the  intervals  |  L . ^ ,  i.e.,  the  amount  of  detail 
in  the  description  of  SN ( R) .  If  the  ensembling  (randomness  in  R(t))  leads 
to  insensitivity  of  detection  probabilities  to  anything  but  the  average 
lengths  of  the{L^},then  a  prediction  of  the  EVA  parameter,  SN,  need  not 
be  very  accurate  in  representing  anything  but  those  average  lengths. 
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SECTION  4 

TRANSMISSION  FLUCTUATION  PROPERTIES 

In  a  previous  investigation  (Ref  4-1)  it  was 
determined  that  the  value  of  the  random-process  (or  any 
other)  approach  to  predicting  system  performance  was  criti¬ 
cally  dependent  on  the  accuracy  of  "input"  data: 

•  mean  transmission  loss  (TL),  at  the  beamformer 
output , 

•  mean  beam  noise  (BN) , 

•  the  temporal  statistics  of  TL  and  BN  fluctua¬ 
tions  (at  a  level  of  detail  consistent  with  the 
amount  of  ensembling,  as  discussed  above). 

Current  acoustic  models  can,  in  many  cases,  predict  the  mean 
values  with  adequate  efficiency  and  accuracy,  given  a  good 
description  of  the  important  environmental  parameters  (e.g., 
sound  speed,  bottom  properties,  noise  sources).  Further¬ 
more,  there  are  models  which  can  estimate  the  temporal 
fluctuations  in  detail  (again,  see  Ref  4-1).  These  models 
are  based  on  the  physical  mechanisms,  and  (at  least  in  the 
case  of  TL)  have  shown  good  agreement  with  measurements. 

Unfortunately,  fluctuation-prediction  routines 
generally  require  substantial  amounts  of  computer  calcula¬ 
tions.  The  reason  is  that  the  models  act  as  "simulators," 
estimating  the  time  series  of  signal  or  noise,  which  in  turn 
are  statistically  analyzed  to  distill  the  important 
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parameters.  In  this  section  and  the  next,  approaches  for 
predicting  certain  key  statistics  of  TL  and  BN  with  minimal 
computational  effort  are  proposed. 

In  summary,  consider  two  methods  for  estimating  TL 
and  BN  properties  for  performance  analyses: 

(a)  Dse  computer  models  to  "simulate"  the  de¬ 
tailed  time  series  of  TL  or  BN.  The  results 
can  be  used  directly  in  sonar  simulations,  or 
else  important  properties  can  be  obtained  from 
time-series  analysis  to  be  used  as  inputs  to 
random-process  models  (which  in  turn  are 
often  used  in  a  Monte  Carlo  sense  to  provide 
replicas  of  the  original  time  series). 

(b)  Apply  existing,  efficient  models  to  estimate 
mean  values  of  TL  and  BN.  Use  consistently 
efficient  calculations  (proposed  below)  to 
predict  important  fluctuation  statistics 
for  use  in  random  process  (or  other)  perfor¬ 
mance  predictions. 

For  any  practical  application,  (b)  should  be  an  obvious 
choice  over  (a).  This  section  and  Section  5  describe  algo¬ 
rithms  for  estimating  TL  and  BN  fluctuation  properties. 

4.1  Transmission  Fluctuations: 

Assumptions  and  Basic  Models 


The  emphasis  here  is  on  passive  sonars,  employing 
narrowband  processing  at  acoustic  frequencies  from  about  50 


to  2000  Hz  or  more  (for  long  range  detection,  the  lower  end 
of  the  frequency  scale  is  most  useful).  Typically  the 
signal  processor  coherently  integrates  (filters)  for  several 
or  more  seconds  and  then  incoherently  averages  the  resultant 
samples  for  several  minutes.  Interest  then  is  in  character¬ 
izing  the  properties  of  the  average  envelope  of  signal  (and 
noise).  [Ref  4-1  discusses  processing  and  detection 
algorithms  for  this  problem] . 

A  stable  CW  source  signal  is  assumed  in  this  study. 
Fluctuations  in  the  received  signal  are  thus  attributed  to 
transmission,  phenomena,  and  a  number  of  mechanisms  are 
possible.  Consider  two  classes:  those  related  to  the 

motion  or  dynamics  of  the  ocean  medium  and  those  related  to 
motion  of  the  source  and  receiver.  There  seems  to  be  some 
concensus  in  the  community  that  when  there  is  significant 
relative  range  rate  between  source  and  receiver,  the  latter 
class  dominates,  and  the  mechanism  is  the  variation  in  the 
received  multipath  interference  field  (see,  e.g. ,  Refer¬ 
ences  4-2  through  4-4).  The  model  considered  here  is  based 
on  this  assumption. 

Following  Reference  4-1,  the  usual  plane-wave  ray 
formulation  for  signal  multipaths  with  fixed  receiver  and  CW 
source  is: 

s(t)  =  £  Anexp[i(wt-knr+4>  )]  (4-1) 

n 

where  An  is  RMS  pressure,  kn  ■  kcos0n,  k  is  wavenumber, 
r  is  range,  w  is  angular  frequency,  0n  is  ray  angle  at  the 
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receiver,  and  <j>  n  is  a  phase  angle  determined  by 
time  and  phase  shifts  along  the  path.  Intensity 
proportional  to 

T 


s( t )s*( t )dt . 


o 


travel 
is  then 


(4-2) 


Current  ray  (and  wave)  models  have  been  shown  able  (at  least 
at  lower  frequencies)  to  predict  quite  accurately  the 
components  of  (4-1),  with  uncertainty  principally  in  the 
relative  phases  (4>n)  of  paths  with  significantly  different 
histories  (see,  e.g.,  References  4-5  and  4-6). 


When  the  source  location  is  fixed,  a  moving  re¬ 
ceiver  simply  samples  the  static  field  of  (4-1)  in  range. 
For  An  and  kn  locally  range  independent  and  r(t)  the 
time-dependent  distances, 

s(t)  Anexp[i(ut-knr(t)+<t>n)]  .  (4-3) 

The  distributions  of  pressure  and  intensity  over  time  are 
thus  determined  directly  from  (4-3)  and  (4-2).  Special 
cases  discussed  in  Ref  4-1  yield  chi-square  (or  gamma)  and 
non-central  chi-square  distributions  with  2  degrees  of 
freedom  for  the  intensity  fluctuations. 


Continuing  with  the  fixed  source  case,  notice  that 
(4-3)  and  (4-2)  give  a  time  series  for  intensity: 

I(t>  *£  kl  +££  Vn  cos  [<Vkn)r(t)  +  *mn1'  <4'4) 

n  m^n 
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1 


I 

I 

where 

^  km  -  kn  =  k(cos0m  -  cos0n). 

A 

The  intensity  fluctuations  are  thus  driven  by  terms  of  form 
AmAnc°s [ ( km-kn ) r ( t )  ] . 

Following  Reference  4-7,  define  a  "spatial  spectrum"  for  I 
as 

^(fyn-kn)  =  AmAn » 

which  identifies  the  dominant  oscillations  of  I  in  range  or 
in  time  (based  on  r(t)).  Just  how  the  component  sinusoids 
combine  depends  on  the  relative  phases  (<)>mn)  and  the  I1 

amplitudes  AmAn. 

The  effect  of  receiver  motion  on  s(t)  can  be 
viewed  in  two  ways,  the  most  obvious  being  that  the  phase  of 
each  path  (n)  becomes  modulated  by  knr(t),  in  general  a 
slowly  varying  function  of  t.  If,  on  the  other  hand,  r(t) 
is  linear  in  t,  say  r(t)  =  vt,  then  the  n-th  component  of 
s(t)  has  constant  phase  and  frequency  given  by: 

,TcosO 

u).  =  (ui-k  V)  =  (1 - - - £).  (4-5) 

k  n  c 

This  is  the  classical  Doppler-shifted  frequency.  Formula 
(4-3)  thus  shows  components  of  pressure  with  individually 
Doppler-shifted  frequencies;  interference  among  the  compo¬ 
nents  is  determined  for  intensity  from  (4-4). 


1 
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While  the  case  of  receiver  motion  is  straight¬ 
forward,  the  source-motion  problem  is  not.  In  fact,  a  very 
large  amount  of  effort  has  been  expended  in  the  investiga¬ 
tion  of  the  difference  bewteen  the  two  problems  and  of 
methods  for  modeling  transmission  loss  fluctuations  induced 
by  the  moving  source  (see,  e.g.  References  4-8  through 
4-12).  Spofford  has,  under  this  contract  effort,  endeavored 
to  determine  the  important  features  of  source-motion  effects 
and  documented  the  results  in  a  separate  report  (Ref  4-13). 
His  findings  are  summarized  next: 

In  the  context  of  (4-1),  the  received  pressure 
for'  a  fixed  source  and  moving  receiver  has  form 

s(t)  Akexp[i(o)kt-wTk)]  , 

k 

while  that  for  a  fixed  receiver  and  moving  source 
is : 


s(t)  VXP  [i«k(t-Tk)). 

k 


U),  =  COS0Q 

K  J  f 

!  +  v  - 5 - 

Here  Tk  is  a  reference  travel  time,  v'  is  radial 
source  speed,  and  0k  is  ray  angle  at  the  source. 

The  moving-source  condition  can  thus  be  modeled  in 
a  range-independent  environment  in  the  same  way  as 
the  moving-receiver  case,  i.e.,  by  sampling  the 
static  field  originating  at  the  receiver  (in  a  re¬ 
ciprocal  sense)  at  the  source's  relative  speed 
(v').  The  errors  thus  incurred  are: 

°  The  Doppler-shifted  frequencies  u)k  and 
Wfc  are  slightly  different  (second  order  in 
v/c) . 

°  The  relative  phase  between  two  paths  is 
distorted,  with  difference 

k*  v*  (Tkcos  ©k  -  TjcosOj). 
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The  precise  locations  of  peaks  and  nulls  are  thus 
perturbed . 

The  first  error  is  known,  and  of  little  conse¬ 
quence.  The  second  is  unimportant  for  most  appli¬ 
cations  . 

In  summary,  the  statistics  of  fluctuations  in 
transmission  loos  caused  by  source/receiver  motion  can  be 
derived  by  generating  the  static  field  from  a  fixed  source 
and  then  sampling  the  field  locally  in  range  at  the  relative 
source/receiver  range  rate.  The  first-order  statistics  and 
spectra  of  the  envelope  fluctuations  can  be  found  from  the 
range-sampled  data  with  error  of  smaller  order  than  v/c. 

4 . 2  Method  for  Estimating  Fluctuation  Statistics 

Under  the  assumptions  of  the  preceding  subsection, 
a  method  for  efficiently  calculating  statistics  of 
transmission-loss  fluctuations  is  presented  below.  Focus  is 
on  the  variance  and  range  periodicities  (spectrum)  of  the 
received  intensity,  employing  an  approach  developed  by 
Spofford  (Ref  4-14).  The  method  accounts  for  certain 
effects  of  processor  averaging  and  beamforming.  A  subse¬ 
quent  section  compares  results  of  the  method  with  statistics 
generated  by  detailed  TL  models  and  discusses  the  sensi¬ 
tivity  of  the  output  statistics  to  environmental  inputs. 

Consider  two  types  of  multipath  interferences 

(a)  Interference  within  an  Order.  An  arrival 
order  is  a  group  of  paths  which  differ  in 
history  by  a  boundary  reflection  near  the 
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source  and/or  the  receiver.  In  the  inter¬ 
esting  case  that  the  source  or  receiver 
is  near  the  boundary,  two  or  more  paths  have 
nearly  identical  travel  times,  and  persistent 
periodicities  in  range  are  observed.  Surface- 
image  interference  (SII)  is  an  example 
of  this  effect. 

(b)  Inter-Order  Interference.  In  this  case, 
paths  of  different  orders  connect  source  to 
receiver,  and  generate  an  interference  field. 
This  is  the  usual  case  for  bottom-bounce 
propagation. 

For  the  first  case  (a),  the  coherent  sum  of  paths  often 
results  in  "slow"  range  variations  which  are  treated  as 
deterministic  and  included  in  the  prediction  of  mean  TL. 
Under  the  assumption  of  constant  sound  speed  between  a 
source  at  depth  zj  ,  a  receiver  at  depth  Z2  and  a 
pressure-release  surface,  the  SII  effect  yields  an  intensity 
of  form 


I  *  4IQsin2(kzisin  0i)sin2(kZ2Sin9i)  (4-6) 

where  IQ  is  the  intensity  for  an  individual  path  and  ©i 
is  angle  at  source  or  receiver  (more  general  results  are 
found  in  Ref  4-15).  The  change  of  0j_  with  range  yields 
periodic  fluctuations  for  I;  the  fluctuation  period  is  large 
when  kzj,  is  small.  When  the  inter-order  case  is  treated 
below,  such  fluctuation  scales  must  be  examined  to  determine 
if  the  slow  variations  are  to  be  viewed  as  part  of  the  mean 
TL  or  not.  The  FACT  model  (Ref  4-16)  includes  SII,  but  not 
inter-order  interference. 
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Turn  now  to  the  more  general  type  of  interference 
(b)  ,  and  recall  equation  (4-4),  rewritten  in  terms  of  Im, 
the  squared  RMS  pressure  or  intensity  for  the  m-th  path: 

E__  1/2 

I  +  7  7  (II)  cos[(k  -k  )r+<p  ]  (4-7) 

m  t-j  v  m  n  n  mn  K"*  • ) 

m  m^n 

The  range-averaged  intensity  is  given  by 


and  the  "relative"  intensity  by 


61 


I 

<!>’ 


(4-8) 


a  measure  oi  the  fluctuations  of  I.  Notice  that  6  I  has  a 
mean  value  of  one  and  variance  (in  range)  of 


These  expected  values  presume  an  average  in  range  sufficient 
to  sample  full  cycles  of  the  cosine  terms  in  (4-7).  When 
this  is  not  the  case  for  certain  paths  (e.g.,  when  slow 
variations  of  type  (a)  occur),  then  their  contribution  to 
the  variance  must  be  excluded  from  (4-9).  On  the  other 
hand,  when  the  signal  processor's  integration  time  is 
greater  than  a  cycle  of  interference  (as  determined  from 
range  rate),  then  corresponding  contributions  should  be 
omitted  in  (4-9).  Such  modifications  are  discussed  further 
below  in  the  context  of  range  periods. 
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Refer  again  to  equation  (4-4)  (or  (4-7)).  As 
noted  earlier,  the  intensity  is  seen  as  the  (phased)  sum  of 
sinusoids  with  range  periods  p  determined  from 


or 


(k  - 
v  m 


V 


mn 


=  2tt 


X 

pmn  cos0  -cos0 
m  n 


(4-10) 


The  relative  "power"  associated  with  this  period  is 
ImIn/<I>^,  and  the  sum  of  all  powers 


(I  E  I  I  )/<I>2  (4-11) 

'  ,  m  n 

m¥n 

is  the  variance  (  cr2)  of  61.  In  other  words,  the  Fourier 
power  spectrum  for  (4-7)  in  the  variable  r  consists  of 
discrete  components  with  periods  P  mn  and  power  ImIn. 

As  usual,  the  variance  of  the  process  is  equal  to  the  sum  of 

the  component  powers  (Parseval’s  Theorem).  If  certain 
components  are  to  be  omitted  (because  their  oscillations  are 
too  slow  or  too  fast,  or  their  angles  are  outside  the 

aperture  of  source  or  receiver),  they  are  removed  from  the 
spectrum  (i.e.,  added  to  the  "d.c."  component)  and  thus  do 
not  contribute  to  the  variance. 

A  computer  routine  has  been  developed  to  perform 
the  calculations  described  above,  including  SII  effects  and 
accounting  for  elimination  of  fluctuation  periods  outside  a 
selected  interval.  A  typical  case  costs  about  $0.50  for 
computer  execution.  The  mechanics  of  the  program  are 

summarized  next. 
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Intervals  of  allowable  range  period  and  angles  at 
source  and  receiver  are  inputs.  A  ray  trace  is  performed 
for  the  local  sound  speed,  so  that  ray  cycle  distance  (X) , 
arrival  angle  (62)*  and  bottom-intersection  angle  (9b) 
can  be  found  as  functions  of  source  angle  (©i).  For  a 
given  range  (R)  ,  the  intensity  (Im( el))  ot  each  arrival 
is  estimated.  Range  periods  from  (4-10)  and  associated 
"powers"  can  then  be  calculated.  When  a  specific  combina¬ 
tion  of  paths  has  period  greater  than  that  allowed  by  the 
preselected  interval,  modifications  are  made  to  ensure  that 
the  power  and  variance  are  correct;  an  example  follows. 

Suppose  paths  1  and  2  have  P12  larger  than 
allowed.  The  component  I1I2  is  then  removed  from  the 
spectrum  and  an  estimate  of  the  coherent  sum  of  the  paths 
made.  In  most  cases  such  paths  are  of  the  same  order,  and 
SI  I  is  important,  so  that  (4-6)  can  be  applied  to  replace 
the  two  (or  more)  paths  by  a  single  new  path  with  intensity 

A 

I^.  Combinations  of  paths  1  and  2  with  those  of  other 
orders  (say,  path  3),  must  likewise  be  combined.  In  the 
computer  code,  paths  are  grouped  by  order  and  the  correction 
of  (4-6)  made  at  the  start  when  appropriate. 

Once  the  spectrum  is  corrected,  the  range  periods 
and  relative  importance  are  directly  available.  At  this 
point,  the  effects  of  the  signal  processor  can  be  incor¬ 
porated — by  eliminating  components  with  range  periods 
smaller  than  allowed.  Finally,  the  variance  is  found  from 
(4-11). 


In  summary,  an  efficient  computer  code  has  been 
developed  to  extract  the  important  statistical  properties  of 
transmission-loss  fluctuations  caused  by  source/receiver 
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motion  in  multipath  conditions.  Examples,  comparison  with 
detailed  TL  calculations,  and  a  sensitivity  analysis  follow 
in  the  next  subsections. 

4.3  Comparison  with  TL  Model  Results 

Predictions  by  the  routine  described  above  were 
compared  with  the  output  of  a  detailed,  fully  coherent 
ray-trace  model.  The  test  case  and  results  are  summarized 
in  this  subsection. 

The  TL  model  employed  is  a  special  version  of 
FACT  (Ref  4-16)  called  "CFACT",  in  which  all  paths  are 
summed  coherently.  Note  at  the  outset  that  the  cost  of  the 
single  computer  run  exceeded  $200,  as  compared  to  the  usual 
FACT  (semi-coherent)  cost  of  about  $1.  and  that  of  the 
fluctuation  routine  of  less  than  $1.  The  details  of  the 
environmental  and  model  data  follow: 

•  The  sound-speed  profile  is  based  on  a  North 
Atlantic  measurement,  and  is  characterized 
by  a  thin  mixed  layer  above  a  weak  thermo- 
cline,  and  pressure  gradient  from  about  800  m 
to  the  ocean  floor  at  3600  m. 

•  The  bottom  is  flat,  with  loss  at  60  Hz  corre¬ 
sponding  to  FNWC  Class  1:  no  loss  to  19° 
grazing  angle,  6  dB  per  bounce  beyond  52°, 
and  a  linear  combination  between. 

•  TL  was  calculated  every  0.04  nm  to  a  maximum 
range  of  130  nm  for  a  source  at  150  feet, 
receiver  at  500  feet,  and  frequency  of  60  Hz. 
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•  To  determine  fluctuation  statistics,  a  2  nm 
intensity  average  of  TL  was  computed  (TL) , 
and  removed  from  TL.  The  resulting  function 
of  range  corresponds  to  the  normalized  intens¬ 
ity 

<5l(R)  =  I(R)/<I>. 

•  <Si  was  sampled  in  range  over  preselected  range 
intervals,  and  the  resulting  sample  sets  used 
to  compute  the  variance  and  range  periods 
(spectrum) . 

The  routine  developed  for  this  study  was  also  used 
to  generate  o  2  and  p  for  the  same  inputs.  Results  of  the 
comparison  follow. 

First,  note  that  since  the  bottom  is  perfectly  re¬ 
flecting  for  grazing  angles  from  0°  to  19°,  there  are 
at  least  two  dominant  bottom-bounce  orders  of  nearly  equal 
intensity  at  most  ranges. 

Figure  4-1  plots  the  dominant  range  periods  versus 
range,  as  precicted  by  the  fluctuation  routine,  while  Table 
4-1  compares  these  predictions  with  CFACT  results;  agreement 
is  good.  Consider  next  that  as  the  range  increases,  the 
number  of  arrival  orders  does  also.  The  result  is  interfer¬ 
ing  arrivals  whose  periods  while  possibly  not  dominant  are 
still  of  some  importance.  To  determine  how  well  the  fluc¬ 
tuation  routine  predicts  these  less  powerful  fluctuation 
periods,  the  range  interval  100-130  nm  is  examined  in 
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Figure  4-1 

Fluctuation-Routine  Prediction  of  Dominant  Range  Periods 


Table  4-1 


COMPARISON  OF  PREDICTIONS  OF  DOMINANT  RANGE  PERIODS 


Range  Interval 

(nm) 

Dominant  p 

for  CFACT  Data 

(nm) 

Dominant  p-Interval 

From  Fluctuation  Routines 

(nm) 

12  -  20 

0.25  - 

0.29 

0.23  -  0.31  at  45% 
of  the  ranges 

25  -  43 

0.25  - 

0.36 

0.30  -  0.44  at  68% 
of  the  ranges 

45  -  70 

0.27  '- 

0.34 

0.21  -  0.34  at  73% 
of  the  ranges 

100  -  130 

0.28  - 

0.32 

0.19  -  0.23  at  80% 
of  the  ranges 
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greater  detail.  Figure  4-2  shows  the  dominating  p's  as 
computed  by  the  routine  at  each  range  step.  The  ■’  s  repre¬ 
sent  the  primary  periods  at  each  range  while  the  x's  are 
secondary  periods  (those  with  power  within  5  dB).  As  noted 
previously,  the  primary  range  period  for  the  CFACT  data  is 
0.29,  about  which  most  of  modeled  results  are  centered. 
However,  the  spectrum  also  shows  fluctuation  periods  whose 
power  is  within  4-5  dB  of  the  dominant  period  at  about 
1.2  nm  and  0.75  nm.  Such  results  are  reflected  in  the 
figure. 


Turn  next  to  the  predictions  of  total  variance. 
Figure  4-3  compares  o2  for  the  fluctuation  routine  (,*•.  ) 
and  the  CFACT  model  (  h— <  )•  The  length  of  the  bar  indi¬ 
cates  ensembling  interval.  Agreement  is  poor  beyond  20  nm 
for  the  following  reason.  The  isospeed  approximation  for 
SI  I  (equation  4-6)  used  in  the  fluctuation  routine  yields 
3-4  arrival  orders  of  about  the  same  intensity,  and  hence  a 
variance  of  about  0.7,  at  many  ranges  beyond  20  nm.  An 
exact  solution,  based  on  the  true  sound-speed  profile  and 
reflected  in  the  CFACT  results,  shows  only  2  orders  (and 
thus  a  variance  of  0.5).  While  the  fluctuation  periods  are 
nearly  unchanged,  their  powers  and  sum  (variance)  are 
significantly  different.  A  routine  incorporating  a  better 
approximation  for  SII  yields  the  x's  in  Figure  4-3,  with 
good  agreement. 

In  summary,  for  a  test  case  with  CFACT  empha¬ 
sizing  interference  of  bottom-bounce  orders,  the  fluctu¬ 
ation-routine  predictions  of  p  and  a2are  quite  satisfactory, 
and  generated  at  a  small  fraction  of  the  cost  of  running 
CFACT. 
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4 . 4  Sensitivity  Analysis 

The  previous  subsections  have  endeavored  to 
show  that  the  fluctuation  model  has  a  sound  basis  in  physics 
and  produces  reasonably  accurate  results  in  comparison  with 
detailed  TL  model  predictions  for  typical  shallow-source/ 
receiver,  low-frequency  scenarios.  Of  interest  now  is  the 
sensitivity  of  model  output  to  the  various  input  parameters. 

There  are  five  basic  inputs  to  the  fluctuation 

model : 


(1)  Frequency 

(2)  Source  depth 

(3)  Receiver  depth 

(4)  Sound  speed  profile 

(5)  Bottom  loss  function 

Although  each  of  these  individually  affects  the 
calculations,  it  is  often  the  case  that  changes  in  several 
inputs  act  in  concert  to  alter  the  output.  To  illustrate 
the  sensivities  of  the  model  and  to  help  develop  a  general 
characterization  of  the  fluctuations,  a  special  case  is 
studied  in  detail  next. 

4.4.1  Special  Case 

Consider  the  scenario: 


(1)  Frequency  ■  25  Hz 

(2)  Source  depth  -  500  ft 

(3)  Receiver  depth  *  60  ft 


(4)  Sound  speed  profile  is  based  on  a  North 
Pacific  measurement.  It  features  a  thin 
mixed  layer,  a  classic  thermocline  to  the 
axis  at  about  1000  m  and  pressure  gradient  to 
the  bottom  at  6000  m. 

(5)  Bottom  loss  is  6  dB  per  bounce  independent 
of  grazing  angle. 

Concentate  on  multi-order  bottom  bounce  propagation.  For 
these  inputs  the  model  produces  predictions  for  fluctuation 
period  and  variance  as  shown  in  Figure  4-4  and  4-5  respec¬ 
tively.  Surface-image  interference  is  an  important  contri¬ 
butor  to  these  results,  and  Figure  4-6  shows  the  isospeed 
approximation  (Equation  4-6)  for  1(6)  under  the  given 
frequency,  source,  receiver  conditions,  but  without  bottom 
loss.  The  angle  at  the  receiver  which  grazes  the  bottom  is 
approximately  10.7°  and  has  a  ray-path  of  35  am.  Shal¬ 
lower  angles  correspond  to  waterborne  paths,  which  are  of 
little  interest  here  because  of  their  long  periods.  On  the 
other  hand,  angles  that  are  too  steep  will  suffer  many 
bottom  bounces  even  at  short  range,  and  be  too  lossy  to 
matter  here.  The  conclusion  is  that  the  dominant  range 
period  (within  a  scaling  factor  of  frequency)  results  from 
Interaction  of  arrivals  with  angles  which  are  small  (above 
10.70)  and  which  correspond  to  peaks  in  Figure  4-6. 

At  the  shorter  ranges  (to  40  or  50  nm)  where  the 
separation  between  arrival  angles  is  greatest,  it  is  likely 
that  their  exists  an  arrival  angle  near  30°  that  has  not 
been  attenuated  too  greatly  by  bottom  interaction.  Also, 
because  of  the  large  arrival  angle  separation  it  is  unlikely 
that  there  are  two  arrivals  near  the  17.5°  peak.  Hence 
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Figure  4-4 

Dominant  Range  Periods  for  Special  Case  (25  Hz) 
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Figure  4-5 

Variances  for  Special  Case  (25  Hz) 
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Figure  4-6 

Image  Interference  Function 


the  "most  likely"  (in  a  sense,  ensembled  over  range)  fluc¬ 
tuation  period  would  be  the  one  corresponding  to  the  inter¬ 
ference  generated  by  angles  near  the  peaks  at  17.5°  and 
30°.  The  period  is 


p  =  — - -  =  0.38  nm, 

cos!7. 5°-cos30° 


as  predicted  (Figure  4-4). 


Now  as  the  range  increases,  the  arrival  angle 
separation  decreases.  Hence  not  only  do  the  arrivals  on  the 
peak  at  30°  suffer  greater  bottom  loss  but  also  it  becomes 
increasingly  likely  that  there  are  two  arrivals  high  enough 
up  on  the  peak  at  17.5°  to  matter.  In  that  case,  the 
period  of  the  fluctuations  is  determined  by  arrivals  at 
approximately  15°  and  20°,  i.e. 

p  =  - h - =  1.25  nm, 

cos  15°-cos  20° 

and  these  periods  increase  in  range  as  the  separation 
decreases.  Again,  this  is  evident  in  Figure  4-4. 


The  other  parameter  of  interest,  variance,  is 
equally  dependent  upon  the  path  intensities  Im  modulated 
by  the  SII  effects  of  Figure  4-6.  Recall  that 
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If  the  Im's  are  all  equal,  then 


ranging  from  0.0  for  one  order  and  0.5  for  two  to  1.0  for 
large  M.  As  noted  above,  at  the  shorter  ranges  there  are 
two  nearly  equal  orders  near  17.5°  and  30°,  local  maxima 
for  Figure  4-6.  Furthermore,  as  the  range  increases  somewhat 
there  are  two  nearly  equal  orders  for  the  peak  at  17.5°. 
But  as  the  angle  separation  decreases  still  further,  addi¬ 
tional  arrivals  have  angles  approaching  17.5°.  Thus,  a  2 
is  approximately  0.5  at  short  range  and  increases  beyond  50 
nm  as  additional  orders  become  important. 

Particularly  interesting  is  the  30-60  nm  range 
interval,  where  the  fluctuation  period  jumps  several  times 
between  values  near  0.4  nm  (caused  by  two  dominant  arrivals 
near  different  local  maxima  of  Figure  4-6)  and  values  at  or 
above  1.2  nm  (caused  by  two  dominant  arrivals  near  the  same 
local  maximum).  As  will  be  discussed  further,  such  jumps 
occur  over  greater  ranges  whenever  coalescing  of  two  arrival 
angles  near  the  same  peak  is  impeded  (say,  for  example,  if 
the  peaks  are  narrower) . 

4.4.2  Observations  About  Model  Sensitivity 

for  the  Special  Case 

With  the  above  discussion  of  the  SII  effect  in 
mind,  the  sensitivity  of  the  fluctuation  model  predictions 
to  variations  in  certain  of  the  parameters  can  be  explained. 
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Frequency  -  As  the  frequency  increases 
(decreases)  the  peaks  on  the  Sll-curve  become 
narrower  (wider)  and  closer  together  (further 
apart).  Thus,  at  a  higher  frequency,  say  55  Hz, 
the  first  peaks  beyond  the  bottom  grazing  angle 
are  at  13.1°  and  18.4°.  Hence  at  short 


ranges 


cosl3 . 1  -cosl8 . 4C 


=  0.60 


(Note:  There  is  also  the  effect  of  scaling 
implicit  in  the  wavelength).  Additionally 
since  the  peaks  are  narrower,  it  takes  longer 
(in  the  range  sense)  for  two  or  more  paths  to 
arrive  near  the  same  peak.  Figures  4-7  and  4-8 
verify  that  the  increases  in  p  and  a2  with 
range  are  more  gradual . 


Source/Receiver  Depth  -  Since  source  and  re¬ 
ceiver  are  interchangeable  here,  suppose  the 
receiver  is  shallower.  In  that  case,  the  SII 
function  consists  of  the  product  of  two  sinu¬ 
soids:  the  one  with  faster  oscillations  is 
controlled  by  the  source;  it  is  modulated  by  a 
sinusoid  with  slower  variations  for  the  shal¬ 
lower  receiver  depth.  As  long  as  the  change  in 
depth  does  not  change  the  bottom-limiting 
property,  model  predictions  of  P  and  o2  are 
relatively  insensitive  to  changes  in  receiver 
depth.  However,  as  the  deeper  source  varies  in 
depth,  the  location  and  width  of  peaks  in  the 
SII  plot  change,  much  the  same  as  for  a  fre¬ 
quency  shift  (the  deeper  the  source,  the  more 
rapid  the  interference). 
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Figure  4-7 

Range  Period  (Same  Case  as  Figure  4-4,  but  at  55  Hz) 
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Figure  4-8 

Variance  (Same  Case  as  Figure  4-5,  but  at  55  Hz) 


Figure  4-9  and  4-10,  compared  with  4-4  and  4-5, 
show  the  small  effect  of  changes  in  receiver  depth  from  60 
feet  to  150  feet,  while  Figure  4-11  and  4-12  show  the 
substantial  changes  in  p  and  a2  when  the  source  depth 
moves  from  500  to  650  feet.  In  the  latter  case,  the  periods 
jump  because  dominant-arrivals  first  appear  at  angles 
corresponding  to  the  same  peaks  in  the  SI  I  curve,  then  at 
different  peaks,  then  back  to  the  same  peak  again. 

To  this  point  all  of  the  examples  have  been 
similar  because  the  angle  at  the  receiver  that  grazes  the 
bottom  corresponds  to  a  local  minimum  in  the  SI  I  plot.  A 
small  change  in  any  of  the  inputs  can  cause  that  angle  to 
occur  near  a  local  maximum  instead.  Then  the  first  peak 
beyond  the  grazing  angle  is  effectively  only  half  as  large. 
The  net  effect  is  much  the  same  as  that  of  increasing  the 
receiver  depth;  namely  it  takes  longer  (in  range)  to  encoun¬ 
ter  two  arrivals  with  angles  corresponding  to  that  first 
peak.  As  an  example,  consider  a  case  identical  to  that  of 
Figure  4-9  and  4-10,  except  that  the  bottom  grazing  angle 
has  increased  from  10.7°  to  12.5°  because  the  bottom 
sound  speed  is  increased  by  30  ft. /sec.  The  resulting 
model  predictions  for  period  and  variance  are  shown  in 
Figures  4-11  and  4-12.  Notice  that  the  periods  tend  to  be 
clustered  about  the  lower  values  of  P. 

Model  sensitivity  to  sound-speed  profile  and 
bottom-loss  function  is  not  quite  as  straightforward  as  that 
for  frequency  and  source/receiver  depth.  A  discussion 
follows. 
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Figure  4-9 

Range  Period  (Same  as  Figure  4-4,  but  Receiver  Depth 


Figure  4-10 

Variance  (Same  as  Figure  4-5,  but  Receiver  Depth 
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Figure  4-11 

(Same  as  Figure  4-4,  but  source  depth  is  650  feet) 


•  Bottom  Loss  Function.  For  bottom  loss  func¬ 
tions  that  are  independent  of  grazing  angle, 
results  similar  to  those  of  the  examples 
discussed  thus  far  can  be  expected.  The  only 
appreciable  difference  is  that  the  higher  order 
arrivals  become  increasingly  important  as 
bottom  loss  decreases;  in  that  case  the  inter¬ 
ference  generated  by  arrivals  near  different 
maxima  in  the  SII  curve  is  significant  to 
longer  ranges. 

Bottom  loss  functions  that  are  dependent  upon 
angle  are  more  complicated.  For  the  important  case  of 
low-loss  bottoms  where  paths  with  grazing  angles  to  15°  or 
20°  suffer  minimal  loss,,  the  computed  values  of  p  tend  to 
vary  more  rapidly  because  there  are  more  arrival  orders  that 
have  not  been  attenuated  by  the  bottom. 

•  Sound  Speed  Profile.  Three  distinct  sound- 
speed  profiles  were  selected  for  comparison: 

-  constant-speed  profile 

the  Pacific  profile  described  earlier  in 
this  subsection,  and 

-  a  pressure  gradient  profile. 
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Although  not  a  general  sampling,  these  three  do  cover  a 
range  of  conditions  and  allow  certain  general  statements  to 
be  made. 


Changes  in  the  profile  should  affect  the  model 
output  in  two  ways: 

•  A  shift  of  all  angles  at  the  receiver  and 
most  importantly  the  angle  that  grazes  the 
bottom, 

•  A  rescaling  of  the  ray  periods  for  the  paths 
with  the  shallow  angles  (the  ones  of  most 
interest) . 

The  second  shift  simply  changes  the  rate  at  which  the 
arrival  angle  separation  decreases  with  range.  The  first 
effect  tends  to  be  more  important:  the  bottom  grazing  angle 
can  be  shifted  with  respect  to  the  maxima  and  minima  of  the 
SII  curve  resulting  in  changes  in  p  and  a 2  as  described 
above.  In  addition,  a  significant  change  in  the  environment 
can  introduce  a  new  maximum  in  the  SII  curve  (on  the  posi¬ 
tive  side  of  the  bottom  grazing  angle)  causing  substantial 
variations  in  p  and  a2.  As  an  example,  compare  Figure 
4-13  and  4-14  for  the  isospeed  case  with  Figure  4-11  and 
4-12  for  the  Pacific  profile. 

A  study  of  the  results  for  the  three  profiles  has 
led  to  a  somewhat  general  way  to  view  the  fluctuation 
models'  predictions  of  p  and  for  bottom-bounce  paths. 
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Consider  four  range  intervals,  labeled  A,  B,  C, 
and  D,  and  the  schematics  for  range-period  and  variance 
given  in  Figures  4-15  and  4-16.  Within  each  range  interval 
there  are  certain  fluctuation  properties. 

A.  Steep-Angle  Region.  Interference  is  dominated 
by  interaction  of  paths  with  significantly 
different  receiver  angles.  The  range  periods, 

P  ,  increases  as  the  first  CZ  is  approached. 
When  SII  is  important,  the  variance  changes 
rapidly  because  of  the  fast  changes  in  arrival 
angles  with  range. 

B.  Region  with  Two  Separated  Orders.  Dominant 
orders  have  angles  near  <f>  i  and  $2  corre¬ 
sponding  to  the  first  two  peaks  of  the  SII 
curve,  so  that  period  is  near 

=  X 

cos<f>^-cos<J>2 

When  kzi  is  large,  there  may  be  several 
important  periods,  corresponding  to  other 
low-angle  peaks  in  the  SII  curve. 

The  variance  cycles  between  that  for  one  order 
(0)  and  that  for  two  (0.5)  as  the  angles  of 
the  orders  approach  and  retreat  from  the 
angles  of  maximum  constructive  interference  in 
the  SII  function.  B's  extent  decreases  as  f 
and  Zi  decrease,  or  as  the  bottom  grazing 
angle  approaches  zero. 


C.  Transition  Region.  Here  the  important  ar¬ 
rivals  are  both  near  <p  i  or  are  separated 
but  near  <j>i  and  $2* 

D .  Region  with  Several  Dominant  Orders  With 
Constructive  Interference.  Beyond  the  second 
CZ  (or  so),  two  or  more  dominant  orders  have 
constructive  SII  near  the  same  angles 

A  lower  bound  on  p,  say  p2,  corresponds  to 
two  paths  separated  by  the  width  of  the 
interference  peak  of  the  SII  curve.  As  range 
increases  more  orders  approach  <Pi  and  each 
other,  so  that  p  increases.  Discontinuities 
are  seen  as  the  result  of  orders  "lost" 
because  of  destructive  interference  when  the 
angles  go  to  zero. 

The  variance  begins  as  that  of  two  orders 
(0.5),  increasing  with  range  as  new  orders 
become  important,  but  discontinuously  drop¬ 
ping  periodically  as  an  order  is  lost. 
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SECTION  5 

NOISE  FLUCTUATION  PROPERTIES 

A  recent  report  (Ref.  5-1)  surveyed  models  of 
beam-noise  fluctuations.  It  concentrated  on  shipping  noise 
and  categorized  the  prediction  routines  as  "Brute  Force" 
(essentially  Monte  Carlo  simulations  over  source  locations, 
etc...)  or  "Analytic"  (direct  evaluation  of  statistical 
expressions  for  input  properties  of  special  form).  The 
search  for  an  efficient  means  to  calculate  noise  fluctuation 
properties,  consistent  with  the  TL  approach  described  above, 
began  with  the  results  of  the  cited  report. 

Brute  force  methods  are  by  definition  detailed 
simulations  and  hence  poor  candidates  for  a  short-cut  method. 
Hence,  focus  was  on  the  Analytic  approaches.  There  are, 
however,  serious  questions  related  to  ensembling,  as  discussed 
in  the  report.  In  particular,  "short-term”  properties  are 
appropriate  for  detection  models,  but  differ  significantly 
from  the  usual  output  of  an  Analytic  model. 

Four  Analytic  models  were  considered: 

USI  (Ref.  5-2) 

Wagner  (Ref.  5-3) 

BTL  (Ref.  5-4) 

BBN  (Ref.  5-5) 

Of  these,  the  USI  model  was  dismissed  for  several  reasons, 
the  most  important  being  the  ensembling  question,  but  also 
because  temporal  statistics  and  TL  details  are  not  accounted 
for.  The  Wagner  approach  does  not  apply  directly  to  beam- 
noise,  although  certain  features  were  found  to  be  useful, 
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as  discussed  below.  Of  the  remaining  two,  the  BTL  model 
was  found  to  be  the  more  advanced,  and  tests  with  data  had 
been  performed.  The  BBN  model  is  similar  to  BTL's.  The 
approach  then  was  to  extract  from  the  Wagner  and  BTL  schemes 
a  method  for  efficiently  estimating  beam-noise  statistics, 
and  then  to  test  the  method  against  Brute-Force  simulations 
from  DSBN  (Ref.  5-6). 

An  important  modifier  for  this  work  are  the  results 
of  Ref.  5-1  using  careful  Brute-Force  calculations  of  omni 
and  beam-noise  fluctuation  properties.  There  it  was  found 
that  a  "chi-square"  model  of  the  distribution  and  an  expo¬ 
nential  autocorrelation  function  were  good  fits  to  the  data 
when  beam-widths  were  not  too  small  and  nearby  ships  not 
overwhelming.  Hence,  the  variance  and  decorrelation  time 
are  the  key  statistics  required.  When  transient  events  are 
important  (e.g.,  single  dominated  ship),  a  "nearest-neighbor" 
approach,  such  as  suggested  in  the  Wagner  model,  as  indicated. 
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